QCD ghost dark energy can not (even roughly) explain the main features of the 

accepted cosmological paradigm 
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We explore the whole phase space of the so called Veneziano/QCD ghost dark energy models where 
the dynamics of the inner trapping horizon is neglected and also the more realistic models where the 
time-dependence of the horizon is taken into consideration. Both cases with and without additional 
non-gravitational interaction between cold dark matter and ghost dark energy are scrutinized. We 
pay special attention to the choice of phase space variables leading to bounded and compact phase 
space so that no critical point of physical interest is missing. It is demonstrated that ghost dark 
energy is not a suitable candidate to explain the presently accepted cosmological paradigm since 
no critical point associated with neither dark matter nor radiation dominance is found in the phase 
space relevant to the model. A transient stage of matter dominance - responsible for the observed 
amount of cosmic structure - is an essential ingredient of the accepted cosmological paradigm. The 
above drawbacks are in addition to the well known problem with classical instability against small 
perturbations of the background density originated from negativity of the sound speed squared. 
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I. INTRODUCTION 

The nature of the present stage of accelerated expan- 
sion of the universe remains as one of the unexplained 
mysteries in physics. Assuming Einstein's general rela- 
tivity (GR) is right leads to a multiplicity of models which 
explain the current inflationary stage of the cosmic evo- 
lution within appropriate data accuracy [HQ- I n spite 
of their relative success none of these models looks like a 
convincing alternative due to several outcomes, such as: 
i) the cosmological constant problems, ii) the coincidence 
problem, iii) stability issues, etc. Ghost dark energy 
(GDE) models belong in this vast gallery pj-flij]. Here 
the cosmological constant arises from the contribution of 
the so called Veneziano/QCD ghost fields 15]. Although 
in flat Minkowski spacetime the QCD ghosts are unphys- 
ical and make no contribution, in curved/time-dependent 
backgrounds the cancellation of their contribution to the 
vacuum energy leave a small energy density p ~ Aq CD H , 
where H is the Hubble parameter and Aqqd ~ lOOMeV 
is the QCD mass scale [H-H [HI, [13] • A very attractive 
feature of the QCD GDE models is linked with the fact 
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that nothing behind the standard model (SM) of parti- 
cles and GR is required to explain the origin of the dark 
energy. There is, however, a serious objection against 
these cosmological models in connection with the sta- 
bility issue. In Ref.[6] the cosmological dynamics of a 
simple model where the GDE energy density is propor- 
tional to the Hubble parameter [5J (we use natural units 
8ttG = Sn/m 2 PL = c = 1) 

Pg de ^aH, a~A 3 QCD , (1) 

was investigated. The authors found that the squared 
sound speed of GDE is negative in the model c 2 — 
dpgde/dpgde < 0, resulting in an instability against small 
perturbations of the background energy density dp oc 
exp (wtlexp (ik ■ r) . The issue was investigated in de- 
tail in [9] where it was found that, due to non-positivity 
of the squared sound speed, both non-interacting and 
interacting GDE models are classically unstable against 
perturbations in flat and non-flat Friedmann-Robertson- 
Walker (FRW) backgrounds. Although several authors 
dismiss the instability problem raised by negativity of 
c 2 by arguing that the Veneziano ghost is not a physi- 
cal propagating degree of freedom and the corresponding 
GDE model does not violate unitarity causality or gauge 
invariance 0, [l(J [Hi, others do not find the argument 
convincing. 

Other aspects of GDE models such as: i) equivalence 
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with kinetic k-essence [10j, ii) connection with Brans- 
Dicke, and f(R) theories, and tachyons [llj, and iii) 
thermodynamical issues 12j, have been also investigated. 
In the latter case the QCD GDE energy density p g d e is 
shown to be related with the radius of the trapping hori- 
zon f T [l2|: p gde = a(l - e)/f T = a(l - e)yfH 2 + k/a 2 , 

where e = fx /2Hfx- If neglect spatial curvature, as we 
do in this paper, the trapping horizon is coincident with 
the Hubble horizon fx = 1/H, and 



Pgde = a(l - e)H, e = -H/2H 2 . 



(2) 



Usually the contribution from time-dependence of the ra- 
dius of the horizon (oc e) is neglected How- 
ever, this is a particular case which, strictly speaking, 
corresponds to dark energy dominance. Hence, this is no 
more than a convenient approximation which can not be 
intended to explain the whole extent of the cosmic his- 
tory. In general, the contribution of e to the dynamics 
has to be taken under consideration as it is done, for in- 
stance, in 14]. In that reference the global dynamical 
behavior of the universe accelerate d by the QCD GDE, 
where p g d e = a(l — e)^jH 2 + k/a 2 , was investigated by 
using the dynamical systems tools. The authors consid- 
ered the general situation when there is additional non- 
gravitational interaction between the QCD GDE and the 
CDM, given by an interaction term - source term in the 
continuity equations - of the form 



Q = 3H(ap gde + bpcdrn), 



(3) 



where the constants a, b, are adjustable free parame- 
ters which, for defmiteness, we shall take to be positive 
quantities. Depending on the values of these free pa- 
rameters two critical points of physical relevance were 
found. Due to the choice of phase space variables in 
14] (SI, = p l /3H 2 ): p = il cdm /n gde , and e, both points 
Pk ■ (j"fe,£fe), k = 1,2, are correlated with CDM/GDE- 
scaling behavior since fj,^ ^ 0, 
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If remove the non-gravitational interaction (a = b = 0) 
one would have, instead, Pi : (0,0), and P2 : (0,1). 
In the non-interacting case the critical points will corre- 
spond to GDE-dominatcd solutions (ft g de = !)• The ad- 
ditional non-gravitational interaction between GDE and 
CDM transforms these critical points into GDE/CDM- 
scaling equilibrium points. In either case there are not 
critical points which can be associated with matter dom- 
inance. Hence, on the basis of this result the QCD GDE 
cosmological model may be ruled out as it is unable to 
explain the formation of structure in our universe. 



Before making any conclusive argument against QCD 
GDE models on the basis of the results of the dynam- 
ical systems study in [l3 |. we note that the study in 
that reference is unsatisfactory due to a problem with 
the choice of the variables of the phase space. Actually, 
if take into account ([2]) we can rewrite the variable p as 
p = p c dm/oi(l — e)H, i. e., p and e do not actually take 
independent values. As a consequence, for instance, the 
critical point P2 may not exist since at e = 1 the variable 
p is undefined in general. Moreover, the variable p above 
is unsatisfactory in yet another regard: it is unbounded 
< p < 00. This means, in particular, that critical 
points associated with CDM dominance - if any - are at 
infinity and may be lost. For that reason, perhaps, the 
authors of the mentioned study apply, additionally, the 
qualitative technique of the so called 'nullcline' to com- 
plement the dynamical systems investigation and they 
found that, for a given region in parameter space, and 
depending on the initial conditions, the end point of the 
cosmic evolution can be dominated by matter. 

In the present paper we aim at a throughout investi- 
gation of the phase space dynamics of QCD GDE cosmo- 
logical models. We consider separately the cases where 
the time-dependence of the radius of the horizon is ne- 
glected (fT]) and when it is taken into consideration (|2|). 
For completeness we study also the impact of additional 
non-gravitational interactions depicted by an interaction 
term of the form © . We pay special attention to a consis- 
tent choice of phase space variables leading, in particular, 
to bounded and compact phase space so that no critical 
point of physical relevance is missing. The results of our 
study will show that, as a matter of fact, since no criti- 
cal point is found which may be correlated with transient 
matter (or radiation) dominance, the more realistic QCD 
GDE models where time-dependence of the horizon's ra- 
dius is taken into account, can not explain the formation 
of structure in our universe. This drawback is in addi- 
tion to the well known problem with classical instability 
against small perturbations of the background density 
originated from negativity of the sound speed squared 



II. REMARKS ON PHASE SPACE ANALYSIS 

Usually the way to test the (theoretical/observational) 
viability of a given cosmological model is through us- 
ing known solutions of the cosmological field equations 
or by seeking for new particular solutions that are phys- 
ically plausible. However, in general, the cosmological 
field equations are very difficult to solve exactly and even 
when an analytic solution can be found it will not be 
unique but just one in a large set of them. This is not 
to talk about stability of given solutions. An alternative 
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way around is to invoke the dynamical systems tools to 
extract useful information about the asymptotic prop- 
erties of the model instead. In this regard knowledge 
of the critical (also equilibrium or fixed) points in the 
phase space corresponding to a given cosmological model 
is a very important information since, independent on the 
initial conditions chosen, the orbits of the correspond- 
ing autonomous system of ordinary differential equations 
(ODE) will always evolve for some time in the neighbour- 
hood of these points. Besides, if the point were a global 
attractor, independent of the initial conditions, the orbits 
will always be attracted towards it either into the past 
or into the future. Going back to the original cosmo- 
logical model, the existence of the critical points can be 
correlated with generic cosmological solutions that might 
really decide the fate and/or the origin of the cosmic evo- 
lution. 

The above interplay between a cosmological model and 
the corresponding phase space is possible due to an ex- 
isting isomorphism between exact solutions of the cos- 
mological field equations and points in the equivalent 
phase space spanned by given variables (x, y, ...). When 
we replace the original field variables H, p c d mi p g de, 
etc., by the phase space variables x = x(H, p c dm, •••)> 
y = y(H, Pcdm, ...), we have to keep in mind that, 
at the same time, we trade the original set of non-linear 
second order differential equations in respect to the cos- 
mological time t (cosmological field equations) , by a set of 
first order ordinary differential equations with respect to 
the variable r = In a: x' — f(x, y, ...), y' — g(x, y, ..), etc. 
The most important feature of the latter autonomous sys- 
tem of ODE is that the functions f(x,y,...), g(x,y,...), 
do not depend explicitly on the parameter r. In 
other words, we are trading the study of the cosmolog- 
ical dynamics of H = H(t), p cdm = Pcdm(t), by the 
study of the flux in r-parameter of the equivalent au- 
tonomous system of ODE. The critical points of this sys- 
tem Pi : (xi, j/i, ...), i. e., the roots of the system of alge- 
braic equations f(x,y,...) — 0, g(x,y,...) = 0, corre- 
spond to solutions of the original system of cosmological 
equations. If a given equilibrium point P a : (x a ,y a , ...) is 
a (global) stable attractor, i. e., small linear perturba- 
tions around P a : x — > x a + £(t), y — > y a + C( T )> ■•■) decay 
with r, 1 then, independent on the initial conditions cho- 
sen x(tq) = x , j/(t ) = j/o, — j every orbit in the phase 



space will approach P a into the future ({r : r > To}), 
i. e., the stable (global) attractor is the end point of 
any orbit in On the contrary, if a given critical point 
P s : (x s ,y s ,...) were unstable, i. e., small perturbations 
around P s uncontrollably grow up with r, then this point 
were a past attractor or, also, the source point of any 
orbit in the phase space. For a third class of critical 
points, the so called "saddle points", depending on the 
initial conditions chosen, orbits in <]/ can approach to this 
point, spend some time around it and then be repelled 
from it to finally approach to the stable attractor (if it 
exists). 2 Suppose we have a typical phase portrait, com- 
posed of a source critical point P s , a saddle point P* , and 
a stable (global) attractor P a . Each one of these points 
corresponds to given solutions of the original cosmologi- 
cal equations, H = H s (z), H — H*(z), and H — H a (z), 
respectively. In the above expressions z is the redshift 
which is related with r: t = — ln(z + I). A also typical 
orbit in the phase space will start at P s for r = — oo, 
then will approach to P* and, after a finite (perhaps suf- 
ficiently long) At, will be repelled by P* to finally be at- 
tracted towards P a . The parallel history in terms of the 
equivalent cosmological dynamics will be the following. 
The expansion starts with a Hubble parameter dynam- 
ics H = H s (z) then, as the Universe expands, the cosmic 
history enters a transient period characterized by the dy- 
namics dictated by H = H*(z). After a perhaps long yet 
finite period Az the cosmic expansion will abandon the 
latter phase to enter into a stage which dynamics obeys 
H = H a (z) lasting for ever. 

We want to underline that in spite of the existing iso- 
morphism between particular solutions of the cosmolog- 
ical field equations of a given model and points of the 
equivalent phase space, most of the exact particular so- 
lutions of the field equations can not be associated with 
critical points. These solutions can be picked up only 
under very specific initial conditions on the phase space 
orbits, i. c., these will be unstable solutions. If, for 
instance, there are not found fixed points which could 
be associated with CDM-dominated solution, this would 
mean that even if this is a particular exact solution of 
the cosmological field equations, it can not describe a 
matter-dominated stage of the cosmic expansion last- 
ing for enough time as to produce the observed amount 
of cosmic structure. This result would dismiss a given 
model since it is not able to explain one of the main fea- 
tures of the accepted cosmological paradigm. 



1 This occurs if the eigenvalues of the Jacobian matrix 

//. fy -\ 
I 9x 9y ■■■ 1 

{■: : J 

all have negative real parts. 



2 Our discussion here is oversimplified since, in general, critical 
points can be of many types, for instance, spiral, etc. Besides, 
there can be found also (un)stable manifolds such as cycles. To 
worsen things there can coexist several local attractors, saddle 
points, etc. 
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III. SIMPLIFIED GDE MODEL WHERE p gde oc H 



In this section we shall investigate the phase space dy- 
namics of the simplified (approximate) GDE model when 
time-dependence of the Hubble horizon's radius is ig- 
nored We start with the simplest situation when 
there is not additional non-gravitational interaction be- 
tween the GDE and the CDM. 



A. Non-interacting model 

For the flat FRW universe filled with GDE, CDM, and 
radiation, the corresponding Friedmann equation is writ- 
ten as 



X — Qgde> y — 



(8) 



After this choice the following autonomous system of 
ODE can be obtained: 



3 — 3a; + y 



-2y 



(9) 



where the tilde denotes derivative with respect to the 
variable r = In a(dr = Hdt). The phase space relevant 
to the present study is given by the following bounded 
triangular region in (x, y)-plane: 



3H 2 = p gde + p cdm + p r , (4) 

where p c & m is the energy density of (pressureless) cold 
dark matter, p r is the energy density of radiation, and, 
Pgde is the GDE energy density assumed to be given by 
Eq. (TTJ) . The energy conservation equations for the differ- 
ent components are: 

Pcdm + 3Hp cdm = 0, p r + AHp r = 0, 

pgde + SHp gde (l + UJg de ) = 0, (5) 

where u> gd e is the GDE EOS parameter. The definition 
for the GDE energy density ([T]), together with Eq.©, 
yield to the following relationship: 



It will be useful to have several quantities written 
in terms of the dimensionless parameter of GDE en- 
ergy density Vt gde = p gde /3H 2 and of the dimension- 
less energy density parameter of the radiation compo- 
nent tt r = p r /3H 2 . In this regard, by taking the time 
derivative of the Friedmann equation and, considering 
equations and (JS|), one is left with 

n/i i \ 3-3flg d e + ft r H 

3(1 + " 3rfe)= 2-n gde = (7) 

Hence, for the GDE state parameter one obtains u> gde = 
— (3 — fi r )/[3(2 — £l gde )], while for the deceleration pa- 
rameter: q = 2 + 3uj gde - 

In order to put the cosmological equations in the form 
of an autonomous system of ODE, we choose appropri- 
ate phase space variables: £l gde and Q r , or, for sake of 
simplicity and compactness of writing, 



* = {(x,y)\0<x< 1, 0<y< 1, 

0<x + y<l}. (10) 

In terms of the above variables, the Friedmann equation 
can be written in the form of the following constraint: 
Qcdm = 1 — x — y, besides, the GDE EOS parameter can 
be written as uj gde = (y — 3)/ [3(2 — x)}. 

Three equilibrium points P Ci : (x Ci , y Ci ) can be found 
in "J, which correspond to different phases of the cosmic 
evolution: 

1. Radiation-dominated phase: 

P r : (0, 1), VL cdm = 0, a- - 1, ^gde = 0. 

This is a decelerating expansion solution (q = 1). The 
eigenvalues of the linearization matrix corresponding to 
this equilibrium point are: Ai = 2, A2 = 1, so that it is a 
unstable critical point (past attractor) in ^. The GDE 
state parameter is u) gde = —1/3. 

2. CDM-dominated phase: 

p m : (0, 0), n cdm = i,si r = 0, n gde = o. 

This phase of the cosmic evolution is characterized also 
by decelerated expansion (q = 1/2). The existence of this 
solution is necessary for the formation of the observed 
amount of structure. The eigenvalues of the correspond- 
ing linearization matrix are: Ai = —1, A2 = 3/2, so that 
it is a saddle critical point in For the GDE state 
parameter we obtain w gde = —1/2. 

3. GDE-dominated, de Sitter phase: 

p dS ■. (1,0), n cdm = o,n r = 0, n gde = 1. 

This late-time phase corresponds to an inflationary so- 
lution (q = —1). The eigenvalues of the linearization 
matrix for P d s are: Ai = —4, A2 = —3, so that the solu- 
tion is a future attractor. This equilibrium point mimics 
cosmological constant behavior since ui gde = — 1. 
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The global structure of the phase space shows that or- 
bits in ^ converge into the r-past towards the radiation- 
dominated stage. In the forward r-direction, depending 
on the initial conditions, these orbits evolve for some time 
in the neighbourhood of the saddle matter-dominated 
point, until they are repelled from this point to, finally, 
converge towards the GDE-dominatcd future de Sitter at- 
tractor. This model correctly describes the fundamental 
stages of the cosmic evolution arising within the presently 
accepted cosmological paradigm: i) a stage of radiation 
domination from which, ii) a period of matter-radiation 
equality and subsequent matter domination emerge, fol- 
lowed by iii) a late-time inflationary stage. 



B. Interacting model 



In the presence of an interaction term Q the equation 
([7]) has to be replaced by the following: 



Pgde 
Hpgde 



Jp = -3(1 +Ugde) 



Hpgde 



(12) 



We use same variables as before © and we obtain the 
following autonomous system of ODE: 



H 



3 - 3x + y - Q/3H 3 



(13) 



Non-gravitational interactions have proven useful to 
overcome the coincidence problem in the context of GR 
models of dark energy. Actually, consideration of non- 
minimal coupling or, alternatively, of additional non- 
gravitational interactions between matter and the dark 
energy has been repeatedly invoked to get late-time 
matter-scaling attractors which allow accelerated expan- 
sion [18l420| . In the case of a quintessence, or a scalar 
field, it is expected to couple explicitly to ordinary mat- 
ter unless some special symmetry prevents or suppresses 
the coupling [ll| [2l|. Besides, it is intuitive to imag- 
ine that, at least in principle, the different components 
of the cosmic mixture may exchange energy-momentum. 
Although there exist stringent constraints on additional 
non-gravitational interactions coming from five-force ex- 
periments, given the yet unknown nature of the cold dark 
matter and of the dark energy, it is possible to hypoth- 
esize a certain exchange of energy-momentum between 
these constituents of the cosmic substratum [22| . 

In the case when additional non-gravitational interac- 
tion between the GDE and the CDM is considered the 
corresponding energy densities no longer satisfy indepen- 
dent conservation laws. They obey, instead 



Pcdm + ^Hpcdra — Q, 

Pgde + 3i?p g de(l + ^gde) = ~Q, 



(11) 



where Q is the interaction term which we take to be given 
by Eq.©. This interaction term is a trivial generaliza- 
tion which comprises several plausible non-gravitational 
couplings studied in the bibliography. For instance, for 
a = one gets Q — 3bHp c d m - Such a coupling arises in 
string theory or after conformal transformations within 
Brans-Dicke theory (23|. If choose a = b one gets 
Q = 3bH(p c d m + Pgde), an interaction term that was 
derived in Ref. |24j while searching for a dynamical reso- 
lution of the coincidence problem. 



where Q/3H 3 = 'iax + 36(1 — x — y). Here the phase 
space is the same as in Eq. ([T0"j) . The magnitudes of rele- 
vance can be written in terms of the phase space variables 
(fl cdm = 1 - x - y): q = [1 - 2x + y - Q/3H 3 }/(2 - x), 
and ujg de = -[x(3 - y) + 2Q/3H 3 ]/[3x(2 - x)]. 
The following equilibrium points are found: 

1. Non-interacting, radiation-dominated phase 

P r :(0,l), Q r = l, Q = 0, 

while uigde is undefined. Since the eigenvalues of the lin- 
earization matrix, Ai = 2, A2 = 1 + 3&, are both positive 
(recall that we are considering here that both a and b 
are positive quantities), this is a unstable critical point, 
i. e., a past attractor which corresponds to decelerating 
expansion solution (q = 1). 

2. CDM-dominated phase 

P m : (0,0), n cdm = l, Q = 9bH 3 . 

As before the GDE-EOS parameter u) g de is undefined in 
this case. The eigenvalues of the linearization matrix 
corresponding to this critical point are: Ai = —1 — 36, 
and A 2 = 3(1 — b)/2, so that this solution is a saddle 
critical point if b < 1, while for b > 1 it is the future 
attractor in Given that q = (1 — 3b)/ 2, for b > 1/3 
the solution is inflationary (q < 0). 

3. CDM/GDE-scaling solution 



p 

Q = 



l-b 



1 + a 
9aH 3 



-A) 



- b 



1 



a 



- b 



, U g de — 



1 



This equilibrium point exists whenever b < 1. In this 
case the energy densities of CDM and of GDE scale as 
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^cdm/Qgde = a /(l — ^)- 3 As long as q — — 1 it is an 
inflationary solution. The eigenvalues of the Jacobian 
matrix at linearisation around P m -sc are: -^l = — 4, and 
A 2 = -[3(1 - 6)(1 + a - &)]/(l + 2a- 6), so that the 
corresponding solution is the future attractor whenever 
it exists. The GDE EOS parameter U! g d e < — 1, hence, in 
this regime the GDE always behaves like phantom mat- 
ter. 

Consideration of an interaction term of the form Q = 
3H(ap c d m + bpgde) modifies the topology of the phase 
space (fTO)) by replacing the GDE-dominated (de Sitter) 
solution in the non- interacting case by a CDM/GDE 
scaling solution. It is curious, however, that at this 
critical point the expansion is de Sitter: H = a(l + 
a — 6)/[3(l — b)] = H . Besides, given the scaling na- 
ture of P m - sc {pcdm/Pgde = a/(l - 6)), and taking into 
account the Friedmann equation, one is left with con- 
stant p cdm = a 2 [a(l + a - 6)]/[3(l - b) 2 } = p cdm ,o, and 
Pgde = " 2 (1 + a - 6)/[3(l - 6)] = p g d e ,o- This seemingly 
strange situation may be explained in the following way: 
usually the energy density of dark matter decays oc a~ 3 
as the universe expands, while the energy density of the 
phantom GDE fluid (uj g d e < —1) should increase with 
the expansion. However, due to the non-gravitational 
interaction Q both fluids exchange energy. In this case 
the flux of energy is from the GDE into the CDM at a 
rate which guarantees that no neat change in the energy 
density of cither fluid occurs as the expansion proceeds. 

We want to notice that, since the scaling (accelerat- 
ingly expanding q = — 1) solution whenever it exists is the 
global future attractor, then, independent on the initial 
conditions chosen, scaling between dark matter and dark 
energy always dictates the late-time cosmological behav- 
ior. In consequence the additional non-gravitational in- 
teraction between the cold dark matter and the ghost 
dark energy does actually smooth out the coincidence 
problem. 

At this point we want to recall that the present model 
of QCD GDE, where p g d e — otH, i. e., where e = 
and time-dependence of the horizon is neglected, is just 
an approximation which is valid only whenever the ghost 
dark energy dominates the cosmic evolution. Hence one 
should not expect that this is a good model to explain 
the whole cosmic history. This means that we should not 
take too much seriously critical points other than the 
GDE-dominated one. As a matter of fact, in the next 
section where the more physically involved model with 
time-dependent horizon is explored, p g d e — cxH(l — e), 
we shall see that only one critical point is found which is 



related with ghost dark energy dominance. 



IV. TIME-DEPENDENT HORIZON 

Within the context of QCD GDE models of acceler- 
ated expansion the contribution from time-dependence 
of the radius of the horizon (oc e) is usually neglected 
j however this assumption is no more than a 
convenient simplification of the model. In general, the 
contribution coming from the dynamics of the horizon 
through the quantity e is to be taken under consideration. 
In the present section we shall be considering how the 
time-dependence of the horizon impacts the asymptotic 
properties of the Veneziano ghost dark energy model. As 
in the precedent section the cases when the CDM and the 
GDE are bounded only gravitationally and when there 
is additional non-gravitational interaction between them 
will be studied separately. 



A. Non-interacting model 

Here the energy density of the Veneziano GDE is given 
by Eq. ([2]) and the dynamics of the horizon is encoded in 
the quantity e = —H/2H 2 . The Friedmann equation ((4]) 
and the continuity equations ([5]) hold true, however, in 
place of ([5]) one now gets 



Pgde 
Hpgde 



Pgde 
Pgde 



= -3(l + LJ gde ) = -— e -2e. (14) 



s Since, according to the observational evidence, at present < 
^cdm/^gde < 1 => < a + b < 1. 



Besides, by taking the derivative of the Friedmann equa- 
tion and considering the definition of e, one gets 4e = 

SiOgdeQgde + 3 + fi r , Or LU g de — [4e — 3 — £l r ]/3£l gd e- 

In order to investigate the phase space dynamics of 
this model it is convenient to introduce the following 
variables: x = a/3H, y = fl r , and e. In terms of 
these variables the Friedmann constraint is written as 
(tt g de = x(l - e)): Qcdm = 1 - x(l - e) - y, while 
q = — 1 + 2e, and ur g d e = (4e — 3 — y)/3x(l — e). 

The following autonomous system of ODE is obtained: 

x = 2xe, y = -Ay (1 - e), 

e' = {[(3 - 2e)x - 4](1 - e) + 1 - y}/x. (15) 

Since < fl g de < 1, and since we are focusing in ex- 
panding solutions only => H < 0, then < e < 1. The 
phase space where to look for physically meaningful crit- 
ical points is defined as the following 3D compact region: 

9 e = {{x,y,e)\0<x<l, < y < 1, 

< e < 1, < 1 - x{l - e) - y < 1}.(16) 
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Only one equilibrium point of (|15[) can be found in 

P d s ■ (1,0, 0), Q ffde = 1, g = — 1, w ff d e = — 1. 

This point corresponds to a de Sitter solution since, e = 
=>■ if = i?o = const. It can be associated with a 
transient stage of the cosmic evolution since, as long as 
the eigenvalues of the linearisation matrix Ai = 2, A 2 — 
—4, A3 = —3, are of opposite sign, Pds is a saddle point 
in this case. 

B. Interacting model 

As before we shall consider an interaction term of the 
form ([3]). The GDE EOS parameter u) g de is given as in 
the former subsection recalling that now £l g d e = x(l — e). 
The continuity equations ([S]) are replaced by while 
(fT4"|) is replaced by 



Pgde 0/-1 1 \ Q 

= -3(1 + (jJgde) - 



-2e. 



Pgde Hp g de 1 — £ 

The following autonomous system of ODE is obtained: 

e' = {[(3 - 2e)x - 4](1 - e) + 1 - y + Q/3H 3 }/x, 
x' = 2xe, y' = — 4y(l — e), (17) 

where Q/3H 3 = 3ax(l - e) + 36[1 - z(l - e) - y]. The 
physical phase space coincides with the one for the non- 
interacting case ty e in Eq.fJTB]). 

The additional interaction between the dark matter 
and the Veneziano ghost dark energy transforms the de 
Sitter critical point Pds into a de Sitter-CDM/GDE scal- 
ing point: 



m—sc • 



1 



1 



-6' 



,0,0 , n, 



dm 



with q = — 1, and uj gd e = — (1 +o — 6)/(l — 6) < — 1. This 
critical point exists if b < 1, a > 0. Notice that, as it 
happened with the approximate model where p g d e = otH 
(subsection IIII B| ). the non-gravitational interaction ([3]) 
produces a very strange situation where the expansion is 
in a de Sitter stage H = H (recall that e = 0), 



3H 2 = 3H% = 



1 + a - b 



1 



-b 



Pgde = 



" Pcdn 



which is fuelled by a mixture of dust (CDM) and a phan- 
tom fluid (GDE with EOS uj gde < -1). Due to the 
non-gravitational interaction there is a transfer of energy 
and momentum from the GDE component into the CDM 
one, resulting in no neat change of the energy densities 

Pgde = Pgdefi and p c dm = Pcdmfi- 



To judge about the stability of the above critical point, 
without loss of generality, we consider a particular case 
when a — b. In this particular case the eigenvalues of 
the Jacobian matrix at linearisation are: Ai = —4, and 
A 2 , 3 = -[1 - 5o ± V25 - 58a + 49a 2 ]/[2(l - a)], which 
show that P m - sc can be either a saddle equilibrium point 
if a < 1/5, or a stable future attractor otherwise (1/5 < 
a < 1). 



V. DISCUSSION AND CONCLUSION 

According to the accepted cosmological paradigm the 
main features of the cosmic history of our universe are: 
i) an early period of inflation, ii) density inhomogeneities 
produced from quantum fluctuations during inflation, iii) 
a flat, critical density, acceleratingly expanding universe, 
iv) energy density budget consisting of roughly 2/3 of 
dark energy (DE) and 1/3 of cold dark matter (CDM), 
and v) matter content: 29% - CDM, 4% - baryons, and 
0.3% - neutrinos. Any successful cosmological model 
should be able to explain the above mentioned features. 
In particular a transient stage of matter-dominance is 
mandatory to explain the observed amount of cosmic 
structure. In this paper we have demonstrated that, as a 
matter of fact, the so called Veneziano (also QCD) ghost 
dark energy model [31417]] can not be a successful cosmo- 
logical model in the sense mentioned above, since there 
are not critical points in the equivalent phase space of 
the model which can be associated with a transient stage 
of matter-dominance. There is only one critical point 
which, depending on whether there are or not additional 
non-gravitational interactions between CDM and GDE, 
can be either a matter-scaling solution depicting de Sit- 
ter inflation (saddle or stable) in the former case, or a 
GDE-dominated de Sitter saddle critical point in the lat- 
ter case. This does not mean that a particular matter- 
dominated solution of the cosmological equations could 
not be picked out under very specific initial conditions. 
What this really means is that in case such a particular 
solution is found it will be very unstable, so that it would 
not last for enough time as to account for the observed 
amount of cosmic structure. 

From the results of the study in section IIIII it can 
be inferred that the simplified QCD GDE model where 
Pgde = otH is suitable from the cosmological point of view 
since, besides the de Sitter point (future attractor) there 
are critical points associated with radiation-dominance 
(past attractor) and cold dark matter-dominance (saddle 
point). Notice, however, that there is no way in which 
the latter points can be retrieved from the more realistic 
model where p g d e = otH{\ — e) in the limit when e — > 0. 
This does not mean that our study is not mathematically 
consistent (it is), but rather that one can not pretend to 
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describe the entire cosmic history with an approximate 
model which is valid only during a particular stage of 
GDE dominance (e = 0). Hence the radiation-dominated 
and CDM-dominated critical points in section IJTll do not 
actually belong in the phase space of the Veneziano ghost 
dark energy model. 

Another interesting feature of the interacting QCD 
GDE phase space is related with the bifurcation phe- 
nomenon. It is linked with sudden changes of the topol- 
ogy of the phase space when jumping from one region in 
the parameter space (a, b) into a neighbouring one. This 
phenomenon is possible since the system of ODE p7|) 
is a two-parameters system of ordinary differential equa- 
tions. There arise certain (critical) curves in the space 
(a, b) which separate regions which are linked with dif- 
ferent, topologically non-equivalent phase spaces. As an 
illustration, notice that the critical point P TO _ SC in sec- 
tion IIVBI exists only if < b < 1. This means that 
in crossing the bifurcation line Bi : {(a, b) — (a, 1)} the 
topology of the phase space suddenly changes: for b < 1 
there exists one critical point, while for b > 1 there is not 



any critical point in Another aspect of the bifurca- 
tion phenomenon in this case is related with the stability 
properties of the critical point P m -sc itself, which depend 
also on the values of the pair (a, b). For more details on 
the bifurcation phenomenon within the frame of QCD 
GDE models see the reference where the issue was 
thoroughly discussed. 

Our conclusion is that there are at least two good rea- 
sons why we should rule out the Veneziano ghost as a 
model for the dark energy: i) it is plagued by classical 
instability against small perturbations of the background 
due to negativity of the sound speed squared, and ii) it 
is unable to explain the formation of structure in our 
universe. While the former drawback above might be 
circumvented through more or less plausible arguments 
d, E3] the latter one is conclusive. 
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